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Summary 


A review of the available literature on the gear housing vibration and noise radiation is 
presented. Analytical and experimental methodologies used for bearing dynamics, housing 
vibration and noise, mounts and suspensions, and the overall geared and housing system 
are discussed. Typical design guidelines as outlined by various investigators will be given. 



A. Introduction [1-7] 


The gearbox vibration and noise are caused by the dynamics of gear tooth meshing 
which can be characterized by the transmission error [1-4]. The transmission error is the 
deviation of gear angular position from its ideal location due to tooth profile and spacing 
error, and elastic deformation of gear teeth and body. Its magnitude is of the order of 
several microns. This action produces gear tooth dynamic forces at mesh frequency, fg m 
(Hz) where 


fgm-Ngf s (1) 

Here N g is the number of gear teeth on the shaft rotating at speed f s (Hz). Several 
harmonics of fgm are also noted in measured data. Additionally one can get side bands at 
fgm - nfs where n is an integer; here f s can represent any shaft frequencies [4-6]. 

These forces excite coupled torsional/axial/transverse vibratory modes of the gear 
shafts and produce lateral and vertical displacements at the support bearing locations. 
Dynamic bearing forces are then generated due to the relative motions across the bearings in 
the radial direction. These in turn cause housing vibration and noise radiation at all mesh 
frequencies. In most cases, the noise radiation from the gear housing is due to flexural or 
bending vibrations of the housing walls [7]. The characteristic of such a wave motion is 
shown in Figure 1 . If the transmissibilities of the mounts and suspensions are high, they 
may serve as paths for the structure-borne noise and vibration from the housing to 
structures attached. These structures will vibrate and/or radiate noise too. This vibration 
and noise generation mechanism for a typical planetary geared system is shown in Figure 2 
[4,8]. The pulsating force form over a one tooth spacing cycle generated at the gear teeth in 
contact for each pair of meshing gears is also shown in this figure. 
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Figure 1 . The flexural or bending wave is also a transverse wave, but in this case the 
motion is perpendicular to both the direction of wave travel and the free surfaces [7] 

There have been numerous efforts since the 1960's to model gearbox dynamics and 
acoustics analytically, empirically, and experimentally. Analytical and experimental 
methodologies have been applied extensively to model the dynamics of geared transmission 
system. Some of these models have included the dynamics of gear housing. However, 
most of the gearbox noise prediction models have been semi-empirical in nature due to the 
complexity of the noise generation mechanism; and there have been many experimental 
programs undertaken to characterize the noise field. The purpose of these studies [7-19] 
have been to predict and control gearbox vibration and noise radiation. Ultimately, the goal 
is to obtain an optimal gearbox design which minimizes its vibration and noise radiation. 

This review presents previous experimental and analytical methodologies used for 
shaft-bearing dynamics, housing dynamics and acoustics, gearbox mounts and 
suspensions, and the overall geared and housing system. These will include a discussion of 
various formulations and assumptions. Typical results and problem areas regarding the 
techniques used will be pointed out. Some of the typical design criteria reported by various 
investigators are also summarized. 
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B. Gear-Shaft-Bearing Dynamics [7-37] 


Gears and shafts vibrations produce bearing reaction forces. These forces are 
responsible for transferring the displacement form of excitations of meshing gears to the 
housing. Knowing the nature of these forces and their transfer paths will allow better 
control and prediction of the gearbox vibration and noise. A detailed review of the gear 
dynamics models has been conducted by Ozguven and Houser [20]. 

Laskin, Orcutt and Shipley [9,10], in 1968, used the Holzer torsional vibration 
models of simple and planetary gearing system to compute gear tooth dynamic forces. A 
segment of this torsional vibration system is shown in Figure 3. Based on this model, two 



I = mass moment of inertia Q = torsional compliance 

T = torque C = damping 

® = angular motion 

Figure 3. General portion of the Holzer torsional system. Subscript n and n+1 indicate 
station number . T91 
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relationships for the angular motion and torque between successive stations and at each 
station were obtained. One equation described the transfer of angular motion and torque 
between stations while the other described the difference between the input torque and the 
output torque at each station. They applied this method, in conjunction with a gear 
excitation model and experimental data, to study the vibration energy paths of the UH-1D 
helicopter transmission. Badgley and Laskin [11,12], in 1970, performed similar 
experimental and analytical studies on the CH-47 helicopter transmission. 

Badgley and Chiang [13-15], in 1972, used a shaft-bearing system dynamics 
approach to obtain the lateral response of a gear support system. Using this approach 
bearing dynamic forces may be obtained from the previously computed gear tooth dynamic 
forces [9-12]. This analysis was performed upon the assumption that transverse vibration 
of the shafts are responsible for transferring the gear tooth dynamic loads to the housing. 
Moreover, the lateral resonance frequencies are within the gear mesh frequency range, that 
is, in the order of a kiloHertz. Finite cylindrical beam elements with rotation and lateral 
degrees of freedom were used to model the system. Nonisotropic linear bearings, and 
uncoupled torsional and lateral motions of the system were assumed. Effects of housing 
flexibility on the gear-shaft dynamics were not included. 

Experimental evaluation was done in parallel with these analytical predictions. 
Qualitative results for gear mesh frequencies, vibration levels, etc. were in good agreement 
with the experimental data. However these methods do not indicate the effectiveness of a 
gearbox design change in terms of vibration and noise reduction in the audible frequency 
range as shown by Stemfeld, Schairer and Spencer [16]. 

Bowes et. al. [17-19], in 1977, reviewed and modified previously constructed 
analytical models by Badgley and Chiang to include the dynamical effects of housing mass, 
stiffness, and damping. Holzer-Myklested technique was used to model the uncoupled 
torsional and flexural vibrations of the geared system with shafts as slender cylindrical 


6 



beams and the gears as lumped masses and inertias. A typical model is shown in Figure 4 
where each shaft segment was treated as uniform torsion-flexure element with distributed 
inertia. The model has N shaft segments with N+l stations. At each station, a rigid body of 



mass m and mass moment of inertia I, and/or an excitation F M may be attached. Due to the 
assumption that the torsional motion is uncoupled from the lateral motion, both the related 
impedance matrices were obtained separately. Using the transfer matrix approach, the 
mechanical impedance of the total geared system was constructed such that [19] 

f, 1 [ ^,,1 U 12 ] [0] 1 rx f I 

F 2 }-= [ Z 21 ] [ Z 22 ] [0] ■| x 2 ( 2 ) 

F 3 J [ (0] T [0] T [z 3J J Kl 



where 


[z .j] = impedance submatrix relating responses in the 

i direction to excitations in the j direction 
[0] = null submatrix 
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x . = responses in the i direction 

F . = excitations in the j direction 

i,j =1,2,3 

with 1 = vertical direction 

2 = horizontal direction 

3 = angular torsional direction 

This resultant impedance matrix was then combined with the housing impedance, which 
will be discussed later, using the component synthesis method. The bearings models were 
nonlinear springs in the two orthogonal directions considering bearing geometry, torque, 
and shaft speed [17-19,21]. The SH-2D helicopter transmission were analyzed [17-19]. 

Salzer, Smith and Welboum [22,23], in 1975 and 1977, simulated a 6 degrees of 
freedom lumped-mass model of an automobile gearbox internal components independent of 
the housing parameter on an analog computer. The system and its simulation block diagram 
are shown in Figure 5 and 6 respectively. An analog model was used. The computed 



Figure 5. Idealization of the automobile gearbox internal components [22] 
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Figure 6. Block diagram of gearbox vibration simulation [22] 

bearing forces were available immediately for audible output through a loudspeaker. The 
results were found to be very similar in character to the experimental data when seen in the 
frequency domain but the magnitudes were not the same. Astridge and Salzer [24], in 
1977, also used partial lumped-mass method to model the vibrations of the Wessex Tail 
Rotor gearbox as illustrated in Figure 7. Thirteen lumped mass locations numbered 1 thru 
13 were selected and each has 6 degrees of freedom. All the shafts and part of the housing 
were modeled as hollow cylindrical beams. However, the stiffness matrix of the complex 
housing section was obtained using the finite element method. Sinusoidal forced response 
analysis indicated very little relative displacement across the bearings since the shafts and 
housing virtually moved together, and the dynamic bearing loads were about 5 % of the 
static loads. 

Recently, Neriya, Bhat and Sankar [25] used the lumped-mass model to include the 
coupled torsional and lateral vibrations of a simple gear-shaft system as illustrated in Figure 
8. At the bearing locations, the simply supported boundary conditions were assumed. 
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A lumped parameter model of dynamometer, motor, shaft stiffnesses, and gears was used 
to obtain a set of second order dynamic equations [9] 

[M]{q}+[C]{q} + [K]{q} = {F} (3) 

where [M] = generalized mass matrix 

[C } = generalized damping matrix 
[K] = generalized stiffness matrix 
{q} = generalized displacement vector 
{F} = generalized force vector 

Using the normal mode analysis method, the dynamic tooth loads were estimated to be 
maximum at the torsional resonances which concluded that coupling between the torsional 
and lateral vibrations did not have significant effect on this behavior. 

The finite element method (FEM) was also used to model the internal components of 
geared transmissions. These models were usually assumed to be uncoupled from the 
housing like most of the above ones by the assumption that the gear-shaft system is much 
more flexible than the housing. Hartman [26] used the finite element method to model the 
transverse-torsional-axial vibration of the 301 HLH/ATC helicopter geared transmission. 
The dynamic tooth forces computed using the approach adopted by Laskin, Orcutt, and 
Shipley [9] were used as inputs in the forced response analysis. He indicated that the finite 
element approach has the advantage of allowing coupling between adjacent shafts across 
the gear meshes by defining gear mesh stiffnesses.The gears were modeled as lumped 
masses and inertias with linear springs between the nodes, shafts were modeled as beams, 
and bearings were modeled as beams and springs. Sciarra et. al. [8,27], Drago [28], and 
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Royal, Drago and Lenski [ 29 ] used similar finite element program to model the CH-47 


helicopter geared transmission as illustrated in Figure 9. 




In addition, the strain energy densities at each mesh frequency were computed to 
identify possible design alterations. Observations of the first 20 mode shapes indicated that 
most of them are primarily coupled bending/torsion modes. The bearing forces computed 
were used to excite the NASTRAN finite element housing model. These loads were phased 
at each mesh frequency due to damping. 

Neriya, Bhat and Sankar [30] specifically studied the effect of coupled torsional- 
transverse vibration of a simple gear shaft system, also shown in Figure 8, now using 41 
degrees of freedom finite element model. A typical beam element with 6 degrees of freedom 
is shown in Figure 10. Nonisotropic bearings elements were assumed by specifying linear 
stiffnesses in two orthogonal directions in the plane of the support bearings. Typical 
stiffness is approximately 10 8 N/m. However, the basis for obtaining the equivalent 
damping coefficient in each mode is not clear. 
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Figure 10. Typical beam element used in rotor dynamics [30] 

Steyer [31], 1987, mentioned that a detail analysis of a geared transmission would 
take up a lot of time and also require some modeling experience. He then suggested an 
impedance analysis of a simple gear-shaft system independent of housing parameters for 
dynamic bearing forces estimation. This was done by assuming a large impedance 
mismatch at the support bearings. First, the excitation at the mating teeth, F mes h> was given 
by the product of the mesh impedance, Zmesh, and the relative velocity between mating 
teeth, ito5 . The mesh impedance was evaluated in terms of the impedances of the shafts 
for the translational and rotational components, and the lateral vibration of the shaft at the 
bearing location was given as [31] 

where Z T1 

(0 


x I = F me S h Z Tf 1/i<0 


(4) 


= mobility of shaft 1 (translational) 


= angular velocity 
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Defining a bearing stiffness, K b ( to ), the bearing force [31] 


F b =x 1 K b (5) 

The final form of the bearing force for identical shaft 1 and 2 was shown to be [31] 

F b = K b 8tZ T (i«>K M ‘ 1 +2Z R " 1 ) + 2r 1 (6) 

where Zt - shaft translational impedance 

Zr = shaft rotational impedance 
Zm = tooth compliance 

The bearing response based on this model is shown in Figure 1 1 . The response was 



Figure 11. Typical force transmissibilitv curve (exact and asymptotic) [31] 
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divided into 5 regions which are also tabulated in Table 1. Each region has its own 
controlling factor, for example, the response in region IV is proportional to the ratio of the 
gear torsional inertia to the sum of inertia and mass, and the response in region HI is 
proportional to the gear torsional inertia 

Table 1. Frequency limits and approximate response for the 5 regions [31] 


Zon« 2fg/K B < 

* «*»* (Jq + JrI/R’Kt 


Frequency Limits 


II 

Kp/IKy ♦ Kp) 

- 

CJ* - 

III 


CJ 1 “ 

IV 

Jq/IR 1 m ♦ J G } 

o 1 “ 

V 


<-> J - 


k t /m 


where 


M - Shaft effective mass 

Ky « Shaft lateral stiffness 

R - Gear pitch radius 

Jq - Gear torsional inertia 

Kr • Shaft torsional stiffness 
- Reaction torsional inertia 


In the previous mathematical models for the vibration energy transfer through the 
bearings, only radial forces were assumed acting through the bearings. Rajab [32] allowed 
radial and moment loads transmitted through the support bearings. The sketch of the ball 
bearing model used is shown in Figure 12. Bearing angular and radial stiffnesses were 
obtained by solving a set of approximate bearing-shaft load-deflection equations using the 
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F m = maximum load 9 = angular deflection 

F r = radial load 5 r = radial deflection 


Mb = moment load r m = pitch radius 



Newton-Raphson iteration method. The solution for the bearing radial force F r , and 
moment Mb were used to define the bearing stiffness elements as [32] 



lbf. / in. 


(7) 



lbf. / rad. 


( 8 ) 
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lbf . in. / in. 


(9) 


K 


0r 





lbf. in. /rad. 


( 10 ) 


These results compared well with the manufacturer data. A typical comparison is shown in 
Figure 13 for the radial deflections for some radial loads. In addition, a review of the 
mathematical models of the bearings is also presented by Rajab [32]. 

• 0*1* MOM INC Macing manumCIIMI 




This model was then used in the building -block system analysis of the shaft-bearing- 
plate model to study the force/motion transmissibility through the support bearings. Related 
experimental studies were performed on a single shaft supported by a flat rectangular plate 


17 



through a radial contact bearing. The plate was clamped at all the edges as illustrated in 
Figure 14. 


Q 



Figure 14. Shaft-bearing-plate setup for bearing transmissibilitv studies [32] 

Taha [33] also analyzed bearing transmissibility using a set of load-deflection 
equations of the shaft-bearing-housing system. The deflection of the housing was taken 
into account when computing the radial and moment loads across the bearings. These 
analysis were only used to study the effect of bearing misalignment on the performances of 
the gearbox such as shaft deflection and bearing life. The Wessex Tail Rotor gearbox was 
analyzed as an example. 

The statistical energy analysis (SEA) method has been used to analyze power flow in 
marine geared transmission from the gears to the housing [7,34]. This SEA approach is 
valid when the modal density is high. A complex system like a gearbox can be divided into 
many subsystems. An energy balance is then performed on the entire system by 
considering energy stored, energy loss to the environment and energy transfer from one 
subsystem to another. The response of each subsystem is computed in terms of the average 
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and standard deviation of the rms response in a frequency band. Lu, Rockwood and 
Warner [34] developed an SEA model of a marine gear-turbine system, using 79 
subsystems and 148 junctions schematically shown in Figure 15, for comparison with the 
finite element method (FEM). The result is shown in Figure 16. 



Figure 15. SEA model of a marine gear-turbine system F341 



Figure 16. Analytical method applicable range [34] 
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The SEA method is obviously preferred in the high frequency range because it is not 
affected by the increase in the number of participating modes as in the FEM. The general 
power flow equation based on the SEA method is given as [34] 

N 

n, + 2n tj -n 2i 

2 

N 

-1,2 n 2 + Sn 2J 

1,3 


T] . = The loss factor of subsystem i 
t| = Coupling loss factor 
E. = Energy stored in the subsystem 

71 !^ = Input power 
<o = Frequency (rad /sec) 

Lyon [7] also used the SEA method to estimate the transfer functions for the energy 
transfer paths in a marine gearbox. He showed that the SEA prediction was better than than 
the lumped-mass model when compared with a 1/4 scale model.These comparison are 
shown in Figure 17. 

In other experimental studies, Ishida, Matsuda and Fukui [35] studied the 
transmission of vibration energy in an automobile gearbox by examining the acceleration 
and noise frequency spectra at various location of the gearbox and its surrounding. A 
schematic of the vibration and noise transmitting paths is shown in Figure 18. It was also 
found that most (95%) of the total gearbox noise came via the structure-borne paths where 


[ -^IN 
where 


^Nl 
" ^N2 


N - 1 

1 'n + I%j 


fE, 


r i 

K. / CO 

1 


in 

E 2 


n 2 / (o 
in 
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E n, 


7t^ / GO 
in 


( 11 ) 
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the fraction Ej, = 95 % was computed as 

e d=V< e s + e a> (12) 

where Ea = output energy density through air-bome path 

Eg = output energy density through solid-borne path 

The output energy densities were computed from the mean noise reductions for the air- 
bome, solid-borne, and total noise. This high structure-borne noise contribution is due to 
the fact that most of the air-bome noise from the meshing gears was reduced by the 
housing. In addition, a free torsional vibration analysis using Holzer's method was also 
performed on this multispeed geared transmission. 

Randall [36,37] suggested examining the vibration data in the cepstrum domain to 
extract certain information on the gearbox vibration which otherwise cannot be obtained 
from the frequency (spectrum) and time domains. The cepstrum is an inverse Fourier 
Transform of the logarithmic power spectrum, or mathematically [36] 

C(t) = [3* {log F(f)}] 2 (13) 

where C(t) = cepstrum 

F(f) = power spectrum of the time signal 

3 *{ } = inverse Fourier Transform 

This cepstrum analysis was reported to allow one to extract periodicity in the spectrum, 
detect increase in sidebands amplitude and spacing which usually implies deterioration of 
geared transmission, analyze spectra of very fine resolution, separate excitation from the 
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vibration transfer path function, etc. Randall [37] used cepstrum analysis to obtain the 
excitation and its transfer path functions from the measured response of a gearbox. This 
could be done because the cepstrum of the measured response being a sum of the excitation 
and its transmission path cepstra. Also the excitation was found to concentrate at higher 
quefrency range as compared to the transmission path function. To show this application, 
consider the spectrum of the measured response [37] 

F(f) = G(f) * H(f) (14) 

where G(f) and H(f)are the excitation and impulse response spectra. Hence, the Fourier 
Transform of logarithmic measured response function in equation (14) is [37] 

3 -1 {logF(f)} = 3 -1 {log G(f)}+3 _1 {logH(f)} (15) 

i.e. the sum of source and impulse response cepstra is the measured response cepstrum. A 
typical cepstrum is shown in Figure 19. The excitation can be seen to dominate at the high 



Figure 19. Measured response cepstrum [37] 


23 







quefrency range. Once the region of low quefrency range where the effects of the impulse 
response is significant is determined, it is possible to curve fit the response to a transfer 
function with known number of poles and zeroes. This can then be subtracted from the 
total cepstrum leaving only the excitation cepstrum. These cepstra can then be transform 
back to the frequency or time domains for diagnostic. 

Lyon [7] performed mode counts on a gear and shaft to study the vibration transfer in 
these structures. He showed that at high frequency, part of the gear-shaft system acts as a 2 
dimensional structure resulting in a higher number of participating modes. For example, 
Figure 20 indicates that the hub of the gear display new circumferential modes, in addition 
to the 1 dimensional shear and bending modes, at frequency above 16 kHz. On the other 
hand, the bending, and inplane longitudinal and shear vibrations of the rim of the gear 
occurs at all frequencies.This occurrence of additional modes result in higher ability in the 
structure to transfer vibration energy. 



Frtqutncy in fcMx 

Figure 20. Mode counts in third-octave bands for rim and hub of a gear [7] 
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C. Housing Dynamics [4,8-11,13,15-19,21,24,26,28,29,33,34,38-48] 


A number of publications [8,9,11,13,16,18,21,26,35,38-45] contain experimental 
data on gear housing vibration due to gear excitation at the mesh frequencies and their 
multiples. Most of these give the transverse acceleration frequency spectra of the housing 
plates. Ishida, Matsuda and Fukui [35], and Lewicki and Coy [44] indicated that higher 
gearbox operating speed implies higher average rms vibration of the housing walls. Also, 
others have realized that the measurement locations significantly affect the measured 
vibration due to change in vibration transfer path function from one point to another. On the 
other hand, housing vibration was found to be quite insensitive to change in geared 
transmission nominal input/output torque. 

Although extensive experimental studies were undertaken, attempts to correlate these 
test results with analytical predictions were limited. One reason may be the complexity of 
the housing geometry involved, for example the CH-47 and UH-1D helicopter 
transmission described in the previous section. To date, modeling of gear housing vibration 
may be grouped as lumped-mass approach, analytical modal analysis, finite element 
method (FEM), and statistical energy analysis (SEA), etc. Some of these methods were 
combined to form a hybrid model and some were aided by other secondary methods in 
order to achieve a simple but reliable dynamic model. 

One of the early effort to model a gear housing as a nonrigid structure where it was 
not coupled to the gear-shaft system was done in 1972 by Badgley and Chiang [13,15] in 
their continuous effort to predict and control helicopter gearbox vibration and noise. They 
applied thin shell theory to characterize the dynamics of finite cylindrical elements of 
variable thickness used in modeling the ring gear housing of the CH-47 and UH-1D 
helicopter transmission. The choice of this element was a natural one for the shape of the 
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gear housings with the ring gear. The CH-47 housing model composed of 3 cylindrical 
shell elements is illustrated in Figure 21. Simply supported conditions were assumed at the 



two edges which allowed only rotation about the circumference. Free and forced vibration 
were performed. In the free vibration analysis, axial and/or circumferential modes were 
found to dominate the behavior as expected. It was noted that although the housing is 
axisymmetric, some modes are not axisymmetric, like the 2nd circumferential mode shape 
in which the amplitude repeats itself twice per revolution shown in Figure 22. An example 
of the first and second axial modes are illustrated in Figure 23. Typical natural frequencies 
of the housing are tabulated in Table 2. Comparison of these natural frequencies with the 
gear mesh frequencies and its multiples indicated that the CH-47 housing would react as a 
forced-response vibration, i.e. no amplification due to resonances, and the UH-1D housing 
would react as a resonant-response vibration. The reason given was that most of the gear 
mesh frequencies for the CH-47 geared transmission were lower than the fundamental 
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Figure 22. Circumferential mode shape (n=2) [13,15] 




First Mode 


Second Mode 


Figure 23. First and second axial mode shape [13,15] 



Table 2. Natural frequencies (Hz) of the CH-47 and UH-1D gear housing [13,15] 


Circumferential 

First Axial Mode 

Second Axial Mode 

Wave Number n 

CH-47 

UH-1D 

CH-47 

UH-1D 

0 

4350 

4380 

13500 

5470 

4 

5220 

4020 

14300 

5370 

6 

6350 

3960 

15450 

5440 

8 

7660 

5800 

16500 

7500 

12 

10950 

9450 

19800 

__ 


natural frequency of the housing, whereas a number of the gear mesh frequencies for the 
UH-1D geared transmission is very close to the first axial with second and fourth 
circumferential modes. 

In the forced vibration analysis, the dynamic tooth loads obtained by Laskin, Orcutt 
and Shipley [9,10] discussed in section B, were expressed as a Fourier series and used as 
the input to this analysis. This exercise could be shown by considering the dynamic tooth 
loads of the form [9] 

F A (0, t) = F A (0) c °s cot (16) 

where F . (9) = circumferential distribution of radial forces 

A 

co = forcing frequency ( rad / sec ) 
t = time 

0 = angular position with respect to gear A (Figure 24 ) 


Figure 24 illustrates the coordinates of the planetary gear system. Expansion of the function 
representing the circumferential distribution of the radial forces, as shown in Figure 25, as 
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a Fourier series led to [9] 


F .(0) = -~ + £ a m cos mO (17) 

Z 1 

where ao and am are the Fourier coefficients. Equation (16) and (17) were used to 
characterize the forcing function due to planet gear A. The dynamic response of the gear 
housing for each Fourier coefficients were then computed, and the form of the response 
function may be written as [9] 

w(0 A , z, t) = b A (0 A , z) cos oat (18) 

Finally, the responses due to planet gear B, C and D were obtained in a similar fashion. 
Using the method of superposition, the total response was constructed by the addition of 
each responses using the appropriate spatial and temporal relationships. Responses of the 
two transmission, shown in Figure 26 and 27, were found to support the prediction that the 
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Figure 27. Normal displacement of the UH-1D housing due to lower planet gear forces [9] 


CH-47 ring gear housing acted as a vibration energy transfer (forced-response) while the 
UH-1D one acted as a noise source (resonant-response). 

With respect to earlier coupled housing and gear-shaft vibrations model, Astridge and 
Salzer [24], in 1977, used the semi lumped-mass approach (the stiffness matrix of the 
complex housing section was obtained using FEM) to model the Wessex Tail Rotor 
gearbox shown in Figure 7. Out of the 13 lumped-mass locations specified as mentioned 
before, each with 6 degrees of freedom, 6 of them are located at the housing structure. 
Although the transmission is quite complex, a simple model was chosen to incorporate the 
dynamics of the gear-shaft and housing into one single model. 

Some experimental methods such as operating motion survey [42,45], and 
experimental modal analysis [32,45,46] were also used to model the vibrational 
characteristics of the gear housing plates and to obtain its system parameters. The 
advantage of using these methods as compared to the purely analytical method is that the 
system matrices are constructed from a response data of a real gearbox, where else the 
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development of an analytical model requires knowledge of the housing dynamic behavior 
and assumptions to simplify modeling procedure. 


The operating motion survey technique involves extraction of the mode shapes and 
natural frequencies by examining the transfer function between 2 points on the gear 
housing. Since this method requires mounting of at least two acceleration measuring 
devices (accelerometers), these devices may alter the system characteristics. Singh, 
Zaremsky, and Houser [45] used this method in addition to structural modal analysis and 
acoustic intensity methods to correlate gear housing plate natural frequencies to their mode 
shapes. The comparison of the second mode shapes using these methods is shown in 
Figure 28a, 28b and 28c. 



Figure 28b. Normalized contours of the 2nd mode (operating motion survey) [45] 
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Figure 28c. Normalized acoustics intensity contours of the 2nd mode [45] 

The experimental modal analysis technique has also been widely used not only in 
dynamic analysis of gearbox but also in many other mechanical systems. Modal analysis 
may be defined as the characterization of the dynamic properties of an elastic structure 
through the identification of its mode shapes and natural frequencies. The general steps 
involve are measurements of force and response signal, determination of frequency 
response function using Fourier Transform, and curve fitting to obtain natural frequencies, 
damping, and transmissibility from one point to another. This method allows one to obtain 
the modes of vibration by avoiding interference from the excitation frequencies. As 
mentioned before, Singh, Zaremsky and Houser [45] used this method to obtain inertance 
transfer function of 75 locations for the housing plate shown in Figure 28. Van Haven, De 
Wachte and Vanhonacke [46] also used the experimental modal analysis technique to 
characterize a gear-motor housing reported to radiate excessive noise. They claimed that the 
fundamental frequency coincided with one of the gear mesh frequency, and by ribbing the 
housing interior shifted the natural frequency away from the excitation frequency. 

Rajab [32] also used experimental modal analysis to model a clamped plate with one 
support bearing on it as shown in Figure 14 of section B. This model together with the 
shaft and bearing models were combined using the building-block system (substructure 
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type) analysis where the total system dynamic matrix is constructed from the individual 
component dynamic matrices. This resultant system matrix equation was used for forced 
response analysis to optimize the bearing location for reduced transverse plate vibration. 


With some advancement in acoustic intensity measurement technique recently, Singh, 
Zaremsky and Houser [45]was able to use this method to perform "in-situ" measurements 
of acoustic intensity veiy close (0.5 in.) to the surface of the vibrating housing plate shown 
in Figure 28. The two-microphone cross-spectrum technique was actually used to obtain 
the housing plate vibration modes which were found to compare well with other methods 
such as modal analysis and operating motion survey as shown in Figure 28. The acoustic 
intensity very near the surface was estimated to be [45] 


where 


with 


I r = (P“,) t 

P! + P 2 
P =— — 



p = sound pressure 
u r = radial velocity 
p = air density 
A = microphone spacing 
{ ) t = time averaged 


(19) 


and where the accuracy depended on the microphone spacing A , and the proximity to the 
radiating surface. 

Bowes, et. al. [17-19], in 1977, as mentioned before in the previous section included 
the effects of housing mass, stiffness, and damping in the gearbox noise and vibration 
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analysis of the SH-2D helicopter transmission. Component synthesis method was used to 
connect the gear-shaft system with the gear housing system. This was done by summing 
the terms in the subsystem impedance matrices which corresponded to the same global 
position. The method used to derived the housing impedance was an incomplete modeling 
technique using modal data and approximate mass matrix. The housing was suspended 
using a low rate stiffness to isolate it from its environment for modal testing. Initially, the 
housing was divided into many elemental masses with its corresponding degree of 
freedom. Bowes, et. al. [17-19] used 44 housing degrees of freedom on the SH-2D 
housing where 20 of which corresponded to the interface degrees of freedom. Then the 
diagonal mass elements were obtained from the elemental masses while the off-diagonal 
elements were estimated. The new modified mass matrix was obtained from the 
approximate matrix by imposing the condition [19] 


where 


{<(>.} [M]|(j>.| = 0 , for i * j 


| <J> . J. = transpose of i-th normal mode 

[M] = mass matrix 

_T<j) = j-th normal mode 


( 20 ) 


In addition, the matrix containing stiffness and damping was computed using [19] 


[K] = [MJ 




[M] 


( 21 ) 
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where 


£2 = i - th natural frequency 
c. = i - th damping coefficient 



The undamped impedance matrix was then obtained from the mass and stiffness matrices 
[19] 

[z] = - co 2 [M] + [K] (22) 

which will be used with the gear-shaft system impedance matrix to analyze the gearbox 
dynamics. 

The finite element method (FEM) was also widely used due the existence of general 
purpose finite element program such as NASTRAN, ISAP-4, SPADAS, ANSYS, etc. In 
most cases, the gear housing was modeled independently from the geared transmission 
with assumed boundary conditions and/or input dynamic bearing/gear forces at the 
interfaces. Kato, Takatsu and Tobe [42] used 480 plate elements on ISAP-4 to obtain the 
vibration modes of a simple gear housing consist of rectangular plates. An example of a 
mode shape computed is shown in Figure 29. 



Figure 29. Vibration mode of a sear housing (0.4m x 0.32m x 0.28m) at 1320 Hz [42] 
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Croker, Lalor and Petyt [47] used isotropic thin flat plate and isoparametric thick flat plate 
elements .shown in Figure 30, on SPADAS to model the vibration of an engine block. 



Substructure method which involve dividing the housing into several parts resulting in 
smaller mass/stifftiess matrices and assembling the global matrices with the assumption that 
each substructure can be adequately represented by only a few modes. Using the properties 
of symmetric and antisymmetric motions, the model size was reduced, but two separate 
analysis were done instead. For example [47], a symmetry about the y-z plane would 
require 

u = e Y = e z = ° (23) 

while an antisymmetric motion about y-z plane would require [47] 

v = w = 0 = 0 (24) 

where u , v , w = displacement in the x, y, z - directions 

0 X , 9 Y 9 z = rotation about the x, y, z - axis 
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A typical correlation between the theoretical and experimental natural frequencies is shown 
in Figure 31. 



Figure 31. Correlation of theoretical and experimental natural frequencies [47] 

In the effort to model the complex CH-47 helicopter transmission housing, Drago, et. 
al [4,8,28,29,48] used the NASTRAN finite element program to develop 3 complex finite 
element models of the CH-47 gear housing parts. The models are for the upper cover, ring 
gear housing, and case as shown in Figure 32. Quadrilateral and triangular homogeneous 
plate elements with membrane and bending capabilities. The 3 sections were analyzed 
separately with simply supported boundary conditions at the interfaces to simulate restraint 
on the boundaries by adjacent sections. Table 3 list some of the natural frequencies of each 
section which are in the vicinity of the planetary gear mesh frequencies. 

Strain energy methods was also used with the above finite element models to calculate 
the strain energy density for each troubled vibration mode. The structural elements with the 
highest strain energy per unit volume were determined as the best choice for structural 
modification. This local alteration of the housing would require minimal weight change for 
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Figure 32. Finite element housing transmission model for CH-47 [8] 

maximum shift in the natural frequency. Areas of high strain energy for modes 3 and 4 are 
shown in Figure 33. 

Finally, the statistical energy approach in characterizing the dynamic behavior of the 
gear housing by statistical means was used by Lu, Rockwood and Warner [34] as 
discussed in detail in the previous section. They summarized that this method is suitable for 
average response determination in the high frequency range. On the other hand, finite 
element method was recommended for estimating the response at the lower frequency range 
due to the detail information available. 
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Table 3. Some of the natural frequencies near the excitation frequencies [48] 


Excitation 

Calculated 

Natural freauencies 

(Hz) 

Frequencies 

UoDer 

Gear 



Cover 

Housing 

Case 

1566 

1518 

— 

1541 


1568 

2334 

1603 

3132 

3069 

2565 

3103 


3133 

3206 

3181 

3606 

3570 

3206 

3588 


3653 

4130 

3664 

4698 

4577 

4130 

4667 


4775 

S — - •* M# ""*■ ■■■“— " 

4770 

4735 


g „ —**» < 





Figure 33. Areas of high strain energy density for modes 3 and 4 [8] 


ORIGINAL PAGE IS 
QE EOOR QUALITY 


D. Noise Radiation [4,7,9-15,17-19,22,23,26,35,38,39,41,42,44,45,49-59] 


Gearbox noise radiation model has been semi-empirical in nature due to the 
complexity of the interactions between a vibrating gearbox structure, such as a gear 
housing and its surrounding fluid. Exact mathematical solution to a sound radiating surface 
in oscillatory motion has been restricted to simple sound sources and highly idealized 
environment, such as a pulsating or oscillating sphere and piston radiator [7]. There were 
many attempts in the past to characterize and correlate gearbox noise frequency spectra with 
the structural vibration and/or excitations spectra using semi-empirical prediction formulas 
and various experimental technique [4,11,13-15,18,19,26,35,38,39,42,44,45,49,50]. 
Most have concluded that the noise prediction is quite complicated and hence an analysis 
requires many assumptions. 

Laskin, Orcutt and Shipley [9,10], in 1968, related the vibration energy in the 
gearbox to noise radiated. They derived a gearbox noise level mathematical expression by 
establishing a semi-empirical relationship between the acoustic energy and gear excitation 
energy. To show this, the total vibration energy Em generated by the gear excitations was 
formulated by Laskin, et. al.[9], as 

E M = ( “ cos (2©t + 0) + 2 cot • sin (0) } (25) 

where 8 Q = excitation amplitude 

F q= force amplitude 
0 = phase angle 

(o = 2it f = frequency of vibration (rad / sec) 
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The first term on the right hand side of equation (25) represented mechanical vibration 
energy, whereas the second term represented dissipated energy through structural 
dampings. The acoustic energy released per cycle E A [9], 


E = a 

A 


5 0 F 0 


(26) 


was then obtained by introducing an energy conversion factor (acoustic efficiency) a for 
the mechanical energy part. By summing all excitations which contributed to the noise level 
at frequency f, the sound power W A expression became [9] 

W A = {af£8 0 F 0 (27) 

Equation (27) was also expressed in sound pressure Lp, dB at distance r by referencing it 
to a standard set of conditions (point source, free-field, atmospheric temperature=68°F and 
pressure=29.5 in-Hg) and introducing geometry and environment factor P . The sound 
pressure level Lp was given as [9] 


L p = 10 log 


Po c o a P f Z 5 o F o 
8 7tr 2 p 0 2 


(28) 


where the reference pressure is p 0 = 2 x 10' 5 Pa and the acoustic impedance is 
p 0 c 0 = 473 kg / m 2 s . It was noted that the accuracy of this formula depended on the 

value of the two factors, i.e. a and P, in equation (28). Badgley and Laskin [11], in 1970, 
rewrote the above sound pressure level expression at third-octave band widths and at full- 
octave band widths by introducing a filter attenuation factor at each band width.These semi- 
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empirical relationships were found to predict poorly when compared to experimental data as 
shown in Figure 34 due to the reason mentioned before, that is the uncertainty in the 
numerical values of the factors involved. However, it was noted from the same figure that 
the character of the noise level is similar to the measured level if the amplitude difference is 
ignored. 



Figure 34. Comparison of empirically prediction gearbox noise and experimental data for 
cruise flight condition (CH-47 helicopter) [11] 

Badgley and Chiang [12,-14], in 1972, estimated the sound power radiated W A by 
the CH-47 ring gear housing using a semi-empirical formula, based on a point source 
assumption and unity radiation efficiency, given below as [12] 

W A = <o 2 w 2 Ap 0 c 0 (29) 

where w = averaged normal displacement of housing 
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co = frequency (rad/sec) 

A = Area 

c 0 = sound velocity 
P o = density of medium 

The computation of the ring gear housing average displacement was obtained from a 
composite cylindrical shell structure model Of the ring gear housing discussed in section C . 
A shortcoming of this formula is no inclusion of the housing geometry. In addition they 
also expressed the equivalent noise level change L^dB as [14] 


L eq (dB) = 20 log 




(30) 


0 N 

where F B and F B are the original and new bearing forces respectively. This changes 
allowed them to evaluate modifications in the geared system design using bearing forces for 
reduced noise level. Similarly, Salzer, Smith and Welboum [22,23] assumed that the 
housing does not change the noise character, and used the bearing force frequency spectra 
to represent the noise in their analysis. 

Section B of this review mentioned Bowes, et. al. [17-19] refined Badgley and 
Chiang models of geared transmission system. In the process, Bowes, et. al. [19] modeled 
the gear housing as a small number of simple, baffled, hemispherical acoustic sources. 
Each source size was estimated as [19] 



( 31 ) 
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where 


s. 

i 


= hemispherical source radius 


R 0 = radius of sphere enclosing transmission housing 
n = number of sources 


Hence based on this assumption, the total sound power radiated (summation of all the 
sources) was computed using [19] 


w A =X 


4 2 2 

S i ( a o). P° C k * 

<o 2 (l+kV 2 ) 


(32) 


where 


( a o) . = absolute amplitude of acceleration at point i 

k = wave number 
Ck = speed of sound 
P o = medium density 


A typical comparison between the theoretical prediction and experimental data is shown in 
Figure 35. Within each frequency band, the prediction closely match the experimental data. 

Ishida, Matsuda and Fukui [35], on the other hand modeled an automobile gear 
housing as a circular piston in an infinite baffle to obtain relationship between the sound 
pressure p and acceleration or velocity of surface vibration. He summarized these 
relationship as [35] 

p « to x when ka<2 or, p «= to x when ka>2 (33) 


where 


x = amplitude of surface vibration (displacement) 



[~~| wwi-aco 

FMOICttO 



k = wave number 
a = radius of the circular piston 
co x = acceleration of vibration 
co x = velocity of vibration 


Hence, the sound pressure level would be either proportional to the acceleration or velocity 
level depending on the area of the vibrating surface. 

The link (radiation efficiency) between the structural vibrations and sound pressure 
level assuming an ideal environment is the most important steps in predicting noise from a 
vibrating surface. It is also very difficult to estimate due to the complexity of the noise 
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generating mechanisms and the fluid-structural interactions, as mentioned earlier. Except in 

very simple cases [51-54], the analytical expression for the radiation efficiency, defined 
below [52], is generally not available. The radiation efficiency 0 ^ is 


a = 
rad 


Po C ° S|v) s t' 


(34) 


where 


= sound power 

S = vibrating surface area 

=mean rms surface velocity (spatially averaged) 
c 0 = sound speed of the medium 
P o = medium density 


Richard [55] realized the elaborate computation and difficulty that one might 
encounter in noise prediction, offered an expression for the A-weighted equivalent sound 
pressure level (db) in terms of structural response, radiation efficiency, damping, machine 
bulkiness etc. given by [55] 

L A, eq (f) = 10 l0 & E escape* 101o g( s - c ) + 10 

- 10 log T) s - 10 log d + B (35) 

I 

) 

where L^eq = A-weighted equivalent sound pressure level 

^escape = total structural energy 

s . c . = fraction of E escape in the frequency band of interest 

A = A-weighted correction 
G rad = radiation efficiency 
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f = frequency (Hz) ORIGINAL ?^ GE B 

5 = damping QE P.OOR QUALI 

d = machine bulkiness 

B =constant 


This formula do not give exact noise level but does indicate the probable factors that might 
explain high noise level in a particular machinery, in this case a gearbox. The contribution 
of each factor to Leq in graphical form is shown in Figure 36. 
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Figure 36. Total noise level in a machinery noise application due to various factors T551 
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Although gearbox noise prediction models has been mostly semi-empirical, there 
exist some numerical methods like the finite element method, finite difference method, and 
boundary element method for noise prediction. These methods are usually difficult to apply 
for complex geometry, and therefore has not been used in gearbox noise analysis. The 
finite element method requires three dimensional acoustic finite element model to 
characterize the noise field exterior to the structure. In addition, there is the problem of the 
termination location for this model which in reality is at infinity for free field condition. 
Therefore, this method is used primarily for closed spaces, and low frequency due to the 
fact that the nodal points spacing must be less than a quarter wavelenght. Finite difference 
method has also similar problem too. Boundary element method [56,57] has been more 
popular because it involves solution to a two dimensional problem of the Helmholtz integral 
equation. It is most suitable for free field sound radiation computation. This method 
requires knowledge of the structural vibration modes which can be obtained using a finite 
element method or an experimental modal analysis. The Helmholtz integral equation is 
given by [56] 


C(y) P(y) = J [P(Q)G'(P,Q) + iz 0 kv(Q)G(P,Q)]dS(Q) (36) 

s 


where Q = surface point 

y = point exterior to the structure 

P = acoustic pressure 

f 2 k if y on surface 
[ 4k if y exterior of surface 

zq = characteristic impedance 

k = wave number 

G = exp (-ikR) / R (Green's function) 
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()'= normal gradient 
v = surface velocity 

The above equation are then reduced to a set of algebraic equation by discretizing the noise 
radiating surface with appropriate elements. These equations will relate normally the 
surface acoustic pressure to the structural surface velocity. 

Few experimental methods such as acoustic intensity method [42,45,52], free field 
measurement technique [12,35,49-51,58,59] and acoustical holography method [42,59] 
were also used widely to characterize gearbox noise level due to the many difficulties 
involved in applying the above methods practically. Free field measurement technique 
requires an anechoic environment whereas the acoustic intensity method allows "in-situ" 
tests. The basis for computing the intensity using this method is given in equation (19) of 
section C. Singh, Zaremsky and Houser [45] used the two microphone "in-situ" acoustic 
intensity method obtain sound intensity very close to the surface of a gear housing plate 
also discussed in section C of this review. Kato [42] performed "in-situ" acoustic intensity 
measurements on gearbox noise in poor acoustical environment. The results indicated that 
certain intensity components intensified by 2 dB (small error) when measurements were 
made near reflecting walls. The explanation given was the diffraction of the sound waves 
occurring. An example of the intensity distribution on the measurement surface around a 
simple 0.4m x 0.32m x 0.28m gear housing is shown in Figure 37. Janssen and De 
Wachter [52] also used intensity method to evaluate the contribution of partial surfaces of a 
housing to the total noise radiated. The information was used to aid in design changes by 
use of blocking mass to reduce noise level. Umezawa and Houjoh [59] has developed an 
acoustical holographic system to show locations of sound sources in machinery. The 
process involved hologram recording, reconstruction of recorded wavefront, and intensity 
distribution calculation. This method was applied to an operating simple gearbox. The 
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results obtained were fundamentally known such as the frequency content, noise source, 
etc. 




Figure 37. Intensity distribution around the simple gearbox obtained using acoustic 
intensity method [42] 
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E. Gearbox Mount System [7,13,51,60-66] 


Basic theory on vibration isolation of simple vibrating system, such as the one degree 
of freedom mass-spring-damper system, has been rigorously treated. The reader is referred 
to references [51,60,61] or other equivalent texts for more information. Here, the mounts 
and suspensions of gearbox will be discussed. As mentioned previously that the gear 
excitations not only caused gear housing vibration and noise radiation but the vibrational 
energy may also be transmitted through the mounts and suspensions to structures attached. 
In addition, there will dynamic interactions between the gearbox mounts and gear housing 
which cannot be ignored. 

One of the earlier attempt to model the helicopter gearbox mounts and suspensions 
was done by Badgley and Chiang [13], in 1972. The model consisted of the gearbox 
mount, isolators, the aircraft structure using combination of mass, linear spring, and linear 
dampers as shown in Figure 38. The isolators was assumed to be massless which resulted 
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Figure 38. Analytical model for the gearbox to airframe isolators [13] 


52 






in only two degrees of freedom system. The vibration source was applied at the gearbox 
mount and was assumed to be oscillatory. Using standard method, the equation of motion 
derived was [13] 


m l X l W X l~ X 2) + C l X l + k o( X l~ X 2) +k l x l = F 0 COS <37) 

m 2 x 2 +C o( X 2- X l) + C 2 x 2 + k o( X 2" X l) +k 2 X 2 = 0 (38) 


where the symbols are defined in Figure 38. Then the force and motion transmissibilities, 
T m and Tf respectively, was obtained as [ 13] 


ryi __ I '-'I 

m ” ,Z 0 + Z 2 


(39) 


T Z 0 Z 2l 

f |z . z _ + z _ z , ■+■ z 


T 


0 1 


0 Z 2l 


(40) 


where 


z 


z 


z 




(mechanical impedance of isolator) 
(mechanical impedance of mount) 
(mechanical impedance of aircraft) 


Based on this simple analysis, Badgley and Chiang [13] concluded that for low motion and 
force transmissibility, the mechanical impedance of the isolator must be small, and the 
mechanical impedances of the local gearbox mount and local aircraft structures must be 
high. In other words, the isolator must be made as soft as possible with low damping, 
while the gearbox mount and aircraft structures must be massive and highly damped. The 
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difficulty in obtaining reliable physical values like mass, spring, and damping was also 
mentioned. Several methods were suggested, that is, numerically compute these physical 
quantities based on the geometry and material, and experimentally extract the impedances. 

Warner and Wright [62], and Andrews [63] investigated various marine gearbox 
mounts and isolators requirements for reduction in the force/motion transmissibilities. 
These studies have resulted in the design of special purpose isolation system. Warner and 
Wright [62] identified the energy source such as the transmission error at mesh frequency 
and unbalance of gears and shafts contributing to the vibration transferred through the 
marine gearbox mounts. The addition of damping at the isolators was recommended to 
damp the rigid body modes which might amplify the unbalance vibration of the gear-shaft 
system, although it may reduce the effectiveness of the isolators. Based on these 
observations, a metallic isolation system, shown in Figure 39 was recommended. 



The performance of these system was not analyzed analytically but was tested 
experimentally. Some of the features of this system included high stiffness, absence of 
creep often occur in elastomeric isolators, compact, etc. Figure 40 indicated results of a 
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free-free test of the vertical isolator. It can be seen to perform as a vibration isolator at a 
very wide frequency range. The ability of the isolator to act as a vibration isolator when 
installed is illustrated in Figure 41. Comparison has been made between the installation of 
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the stiff steel connectors and the metallic isolation system. Reduction in the structure-borne 
vibration is observed for the case with the metallic isolators installed. 

Andrews [63] utilized the one degree of freedom system isolation concept as a basis 
for the gearbox mount dynamic model. The gear housing and subbase for the entire system 
were assumed to be rigid. Only vertical motion was allowed in the mount model and the 
journal bearing was modeled as linear stiffness. Modal analysis of this system using the 
model described here indicated that the first two modes were shaft deflection type, and the 
third and fourth modes were associated with the vertical motion of the mounts. This 
analysis led to the design of an isolation system shown on Figure 42. The two side 
rectangular blocks were attached to the gear housing while the middle on to the subbase. 
Two isolators, one on each side, were required to mount the marine gearbox. Application 
of this design led to lower gear housing vibration and equality of bearing loads. 


ISOLATION MOUNT 



W ¥I(W 



ttac vkv 

Figure 42. Gearbox isolation system [63] 


Snowdon [61] also discussed in detail characteristics of damped discrete and 
continuous vibration isolators such as elastomeric isolators, combination of spring-damper 
system isolators, and rods-beams system isolators. The examples were not specifically on 
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gearbox application but more towards general machinery application. Lunden and Kamph 
[64] investigated numerically and experimentally the vibration characteristics of a 
lightweight skeletal machine foundation (grillage) as continuous system isolator. They 
concluded that by applying "blocking mass" and damping (discrete and distributed) on the 
system will result in reduction of grillage vibration over a broad frequency interval, and 
hence lower transmissibility through the grillage system. The damped second order 
Rayleigh-Timonshenko beam has been used in the numerical studies. 

Granhall and Kihlman [65], in 1980, expressed the need for knowing structure-borne 
sound sources data of a machinery in order to aid in design of mounts and isolators and 
noise predictions. For this reason, they analyzed a one dimensional vibration isolator 
system using the mechanical impedances in an analog circuit, and formulated an equation 
for estimating insertion loss of an isolator from measured impedance data. The insertion 
loss IL [65] is given by 


IL =20 log 


z f z i +z f z m +z m z j 

Z i( Z f + Z m) . 


(41) 


where z m , zf, and Zj are the internal, foundation, and isolator impedances respectively. If 
one assumed that the foundation is very rigid, equation (41) may then be written as [65] 

IL = 20 log [l + -t 2 -] (42) 


Comparison of the insertion loss predicted by equation (42) with measured insertion loss 
data, and the insertion loss of a mass-spring-damper system model is shown in Figure 43 
[65]. The graphs indicate that equation (42) predicts the measured data better than the one 
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predicted by a spring-mass-damper model. However, these results are not found to be true 
at high frequencies where both models are inapplicable. 



FREQUENCY, Hz (b) 


Figure 43. Insertion loss in 1/3 octave bands for a fan unit ( solid lines = measured data, 
dashed lines = equation (42). and dotted lines = mass-spring system! T651 


Unruh [66] developed a finite element dynamic model of an aircraft engine mount to 
be coupled with the rigid engine model, frequency dependent stiffness model of the 
isolators, and experimentally obtained fuselage and interior response model. The purpose 
was to study the effect of isolators and mounts on the structure-borne noise transmission. 
The vibration isolator modeled as frequency dependent radial k R and axial k A springs in 
local coordinate was given as [66] 


k R = k R *(o>) [1 + iT](co)] 

(43) 

k A = k A* (a) H 1 + h 1(®)] 

(44) 
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where 


k * = radial spring modulus amplitude 
k A * = axial spring modulus amplitude 

T| = material loss factor 


The finite element model of the mount system, illustrated in Figure 44, consisted of 70 
elastic beam elements with 201 degrees of freedom. Using modal synthesis method, as 



described in previous section, the number of degree of freedom was reduced to 51 elastic 
and 6 rigid body degrees of freedom. For each of the subsystem listed above, the standard 
second order differential governing equation was derived. Then by proper choice of the 
independent degrees of freedom, each components were coupled together by the 
summation of interface forces and was set to zero to obtain an empirical relation between 
the structure-borne noise at various position in the aircraft interior and the chosen degrees 
of freedom. 

Lyon [7] also performed similar analysis on a marine gearbox system schematically 
shown in Figure 45. This method involved modeling of the gearbox mount system in detail 
using combinations of simple beam, spring, damper and mass elements. The input and 
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Figure 45. Marine gearbox mounted on a foundation which sits on isolators [7] 

transfer impedances of all the elements were assembled into a complete system according to 
the numbered nodes while setting the total force at each junction equal to zero or to the 
externally applied force. The impedance of these simple elements can be derived easily. 
Figure 46 illustrated the model of a reduction gearbox mount system. A set of mass 



Figure 46. Model of reduction gear mounting system [7] 
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elements were used to represent the gears, and the case rail was used to model the 
foundation structure also shown in Figure 45. The two system rested on a set of spring- 
damper isolator mounts. All these were then supported by a massive beam structure 
(subbase) which in turn sat on the hull elements modeled as sets of springs and dampers. 
The cross section of the case rail and subbase are shown in Figure 47. It was also noted 
that this technique is very similar to the finite element method except here the transfer 
function used to define the elements are functions of frequency. 
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Figure 47. Cross section of the case rail and subbase [7] 

The purpose of the above studies on gearbox mounts and suspensions was to obtain 
parametric design values that will lead to, lower force/motion transmissibility. In most 
gearbox noise and vibration analysis, the mounting system was not taken into account due 
to the complexity of the gearbox mounts especially in aircraft where the structures are 
geometrically complex. It is also coupled dynamically to the gearbox and fuselage. 
However, the inclusion of the mounting system into the dynamic model is a must to obtain 
a noise and vibration prediction models that truly represent the operating conditions of a 
gearbox. 
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F. Overall Gearbox Dynamics [16-18,67] 


Noise and vibration prediction and control ideally require an analytical model of the 
entire gearbox system, its attachments, and other structures connected (i.e. aircraft 
fuselage, subbase, foundation, etc.in an aircraft application). This is due to the fact that the 
dynamics of each components which serve as vibrational energy paths may have significant 
effects on the overall system dynamics. For example, the low to high discrete frequencies 
excitation generated by the meshing gears in an aircraft are transmitted to the airframe 
through various structural paths such as the shafts, bearings, housing, mounts, and other 
attachment points. Discussions in the previous sections of this review have indicated that 
the dynamics of these structural paths are important to the understanding of the overall 
dynamics. There is nothing in the literature that offers a rigorous treatment on the overall 
gearbox dynamics which includes dynamic interactions between the gear-shaft system, 
support bearings, gear housing, gearbox mounts and suspensions system, and noise 
radiation. Although, there is a need for such model, many difficulties like allowable model 
size for computer implementation, complexity of the noise generation mechanism, dynamic 
coupling between gearbox components, etc. hinder the development of an ideal model. 
Hence, in most cases one or more components are modeled in detail, and the other parts are 
modeled with only few degrees of freedom or assumed uncoupled from the rest of the 
gearbox. These assumptions often limit the applicability of the analysis to specific type of 
gearbox model such as those discussed in the previous sections. 

Berman [67] pointed out the difficulties involved in having a complete dynamical 
model of the gearbox and fuselage. Some of the problems addressed here are: 

1 . Cost involved with the assessment of parametric variations 

2. Inadequacy of finite element models in the acoustic frequency range 
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3. High frequency content of the excitations which often excite many modes of the gear- 
shaft and gear housing system, and so large number of degrees of freedom are needed 

4. Complexity of gearbox geometry that is difficult to incorporate, especially in modeling 
techniques other than finite element methods 

5. Difficulty in modeling interface components analytically 

6. Problems associated with combining various gearbox component models to form a 
complete dynamical model 

In view of these problems, Berman [67] presented a methodology to be used in the 
complex gearbox system. It include independent component representation, improvement 
and development of analytical model using test data, coordinates reduction in the frequency 
domain, component coupling, and implementation on a computer. In component 
modeling, each components may be modeled separately using whatever appropriate 
techniques available, for example finite element model for the gear housing , experimentally 
obtained impedance matrix to represent the fuselage dynamics, etc. By doing so, each 
model may be modified without changing the other components.which allows evaluation of 
design modification to be done easily These models are used with reduced degrees of 
freedom to synthesize the complete gearbox model in the frequency domain of the form 
[67] 


( [K] - co 2 [M] - ito[C] ) X(o>) = F(co) (45) 

where [K], [M], [C], are the stiffness, mass, and damping matrices respectively with 
F( to ) as the excitation vector. The reduced component model retains only the interfaces 

and points of applied force degrees of freedom which usually led to significant reduction in 
the degrees of freedom. This step of reduction in the degrees of freedom can be shown by 
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considering a component impedance matrix reordered such that the retained degrees of 
freedom are in the submatrix zy [67] 



With some manipulation, the reduced impedance Z R becomes [67] 

Z R (CD) = Z l- Z 2 Z 4“ lz 2 T 


(46) 


(47) 


Finally component coupling can be done by the summation of all the relevant degrees of 
freedom in each components, for example if an interface displacement vector xj is related to 
the displacement vector X of the complete system by the expression [67] 

x . = T . X (48) 

i i 

where Tj is the transformation matrix, then the impedance matrix Z( to ) of the total system 
would be [67] 


z (cD)=X T i Tz i T i (49) 

i 

A summary of this method is shown in Figure 48. This method was used by Bowes et. al. 
[17-19], also discussed previously, to model the SH-2D helicopter transmission. The 
analysis was not entirely analytical, for example the gear housing impedance was derived 
experimentally due to the problems discussed here and elsewhere. Also there were many 
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assumptions like the simple radiation model which did not include environmental effect and 
housing geometry, and omitted the effects due to gearbox mounts and suspensions. In 
other gearbox analysis, similar problems also arise. One major difficulty is to be able to 
model the interface components, such as support bearings, gearbox mounts and 
suspensions, attachments with simple models yet detail enough to include significant 
dynamical effects on the entire gearbox system. 
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(i. Design Guidelines [4,5,8,13,14,16-19,29,31,49,52,54,68-74] 


Badgley [14] reported that gear mesh excitations are present even in very high quality 
gears which can be amplified by the resonances in the gear-shaft and gear housing systems. 
Hence, the vibration and noise sources control alone is not sufficient. In order to effectively 
control gearbox vibration and noise, design changes in the force/motion transfer paths, i.e. 
gear body, shaft, support bearing, gear housing, gearbox mount and suspension, and 
connected structures are inevitable. Also, it is worth mentioning that design modifications 
in gearbox are very dependable on the gearbox environment, and its application such as 
helicopter or industrial transmission. 

Some design guidelines for noise and vibration control of the gearboxes have been 
developed in the past. Lack of comprehensive design criteria and proper evaluation 
techniques have resulted in a number of conflicting requirements as suggested in the 
literature. This section presents some relevant design criteria for various components of a 
gearbox other than the gears for reduction in vibration and noise. 

G.l. Gear Support System 

If the shafts are found to have high amplitude of vibration, stiffening parts of the 
shafts may reduce the amplitude of vibration especially at the support bearing locations 
where the force are transmitted to the housing [13,14,17-19,31,68]. This can be done by 
adding mass around the shafts without increasing the mass center offset, or using materials 
with high modulus of elasticity - essentially by changing the natural frequencies of the gear- 
shaft system [29]. It is desirable to have the excitation frequencies away from any natural 
frequencies as it would be in any design. An example of successful implementation of shaft 
modification by the addition of mass is shown in Figure 49 where the amplitude of 
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Figure 49. CH-47 transmission shaft vibration amplitudes for nominal and modified 
configurations [13] 


vibration is reduced significantly. Route [69] suggested that when designing a geared 
transmission system, the highest degree of stiffness permitted by size and weight 
limitations should be specified. 

An alternate method to minimize the force/motion transfer to the housing is to locate 
the support bearings at the node points on the shafts [4,49,68], and/or support the bearings 
using stiff frame [29]. By increasing the bearing stiffness with proper choice of bearing 
type will increase the natural frequencies of the system which may be useful [14,29,31]. 
Drago [4] noted that gearbox noise level usually decrease with increasing preloads. 
However, adverse effect may occur in other areas of mechanical design. Figure 50 
indicates effect of the shaft support bearings system on the overall noise level. Sleeve 
bearing are recommended for use as support bearings in gearbox. Although tests have 
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indicated that the bearing quality in terms of noise reduction is as shown in Figure 50, care 
must be taken when using such a guideline due to the fact that the performance of these 
bearings depended on the other gearbox components too. That is, the type of bearing 
installed will have different effect on the overall gearbox system dynamics by altering the 
natural frequencies and vibrational energy paths. 


Some bearings esn reduce drive noise 



). Effect of various bearings on the overall geared transmission noise [4] 


Filling of hollow shafts with damping materials is also helpful in reducing the 
dynamic response of the gear-shaft system when resonance condition exist [14,16,70]. 
Stemfeld, Schairer and Spencer [16], and Drago [4] tested the effect of damping 
(elastomeric material) applied to gear body on the overall vibration and noise level. The test 
results indicated some vibration reduction occurs but not enough to be used alone in 
design. Hence, it may be used as a supplement to other design changes. Other than the use 
of damping to absorb vibration, use of vibration absorber has also been suggested to 
attenuate vibration in the gearbox. The idea of a vibration absorber is that when the 
absorber is properly tuned, the structure attached stops moving at a particular excitation 
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frequency. This concept is illustrated in Figure 51 where mode 1 shows the in phase 
vibration of the absorber and structure at some frequency and mode 2 shows the out of 
phase vibration with respect to each other at a higher frequency. Hence, somewhere in 
between at the tuned frequency, the structure will stop moving. Again tests performed on 
the absorbers indicated that only some reduction in vibration is observed but not 
significantly to be used alone in design. This is due to the fact that the vibration absorber 
works only at a particular excitation frequency which is usually varying over a small range. 
Moreover, there are mesh frequency sidebands which are not attenuated since the absorber 
is timed to the mesh frequency only. 
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Figure 51. Concept of dynamic vibration absorber [16] 
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G.2. Gear Housing and Gearbox Mounts 


The gear housing is the major noise radiator and also serve as a path for the bearing 
excitations to the gearbox mounts. Selective stiffening parts of the housing will reduce its 
vibration amplitude and increase system natural frequencies [4,8,29,68,69,71]. The 
method used in selecting probable locations for modification in stiffness and mass is 
discussed in the gear housing dynamics section. The basic idea is to perform a finite 
element model of the gear housing to identify its natural modes. Then for each mode, the 
strain energy density is computed and regions with the highest energy density will be 
selected for this process [4,29] as shown in Figure 33 (section C). This approach allows 
minimal change in mass and stiffness of the entire gearbox to achieve an amount of increase 
in natural frequencies. 

Over higher frequencies where the radiation efficiency is almost unity, addition of 
damping through viscoelastic material, and restraint on the gear housing will reduce the 
mean rms transverse velocity of the housing plate and hence the sound pressure level too 
[5,29,49,52,72,73]. Effect of various reinforcements added to the a ring gear housing is 
illustrated in Figure 52. It shows higher reduction in the response for center and ends 
reinforcement applied together than when applied separately. However, this may not be 
always possible due to the weight penalty imposed. Addition of mass on the application 
point of external force, also known as blocking mass method, has shown to reduce noise 
intensity level of a gearbox as seen from Figure 53. 

Some other undesirable gear housing geometries are large flat areas and gently curved 
surfaces because they usually vibrate freely and are good noise radiator. One way to reduce 
these effects are to decouple the areas by slotting the housing, add dampers, and thickening 
the housing [4,5,68]. If weight is not a constraint in the design, the use of cast iron which 
has good sound absorbing properties is recommended [4]. In terms of structure -borne 
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Figure 52. Comparison of ring gear housing vibration amplitude with various 
reinforcement configurations [13] 
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Figure 53. Intensity spectrum of an in-line gear unit, before and after application of 
blocking mass [52] 

paths function, it is better to always supply a rigid load paths between the support bearing 
locations on the housing and attachments points for the gearbox mounts to reduce housing 
vibration. Isolators are used to provide resilient support for the gearbox and to reduce 
force/motion transmissibility through the mounts [54,72,73,74]. This is most useful in 
marine and industrial type application since a massive foundation can be provided. When 
designing a mount-isolator system for reduction in force/motion transmissibility, it is 
desirable to have high mounts and foundation impedances, and low isolator impedance 
[ 13 ]. 
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H. Areas of Further Research 


This review indicates that gearbox dynamics and acoustics pose a major problem in 
the development and implementation of gearbox system technology. The literature confirms 
this as Mark [1,3], Badgley [11-15,49], Bowes [17-19], Drago [28,29,48,68], Ishida 
[35], and others have concluded that gearbox noise and vibration levels in aircraft, 
automobile, etc. are often higher than the allowable limits with respect to human comfort, 
and machinery failure and life. These problems become more acute at high gearbox 
operating speed which give rise to excitation frequencies in the order of several kiloHertz as 
in the aircraft gearbox application. Although many attempts were made to characterize the 
dynamics of gearbox system components, currently no comprehensive design criteria exist. 
Moreover, the literature contains conflicting reports concerning relevant design guidelines. 
All these are mainly due to a lack of the complete understanding of the vibration and noise 
generating mechanisms of a gearbox system. Hence, further research on gearbox dynamics 
and acoustics is required. 

A major portion of the gear excitation energy is transmitted through structure-borne 
paths. However, it is difficult to represent the force/motion transfer through the gearbox 
system analytically and obtain reasonable prediction of the gearbox components vibrational 
level . It would be useful to be able to characterize the transmissibilities, and to identify the 
paths quantitatively. 

Also, in order to successfully derive the force/motion transfer model, the dynamics of 
each gearbox components must be known. The bearing subsystem is yet to be modeled 
with success experimentally or analytically. In addition, the bearing interface models are 
sometimes difficult to characterize due to its compliance, and the requirement of matching 
boundary conditions and continuity at the interface. This will lead to better gear-shaft- 
bearing-housing models. 
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Another major area which is not well understood is the effect of mounts and 
suspensions on the force/motion transmissibility and gear-shaft-bearing-housing dynamics. 
In most gearbox application especially in aircraft, the gearbox are mounted resiliently onto 
the airframe which is usually light and flexible. Here the vibration is found to be 
excessive. 

The prediction of the noise radiated by the housing and other structures attached will 
remains a major challenge. This require a model that can relate the structure vibrational level 
to the sound power radiated. To summarize, the areas related to gearbox dynamics and 
acoustics which are not well understood are: 

1 . Bearing dynamics and interface modeling 

2 . Force/motion transmissibility study including an evaluation of the energy paths 

3 . Gearbox mount and suspension dynamics and their effects on the overall dynamics and 
acoustics 

4 . Noise radiation prediction from housing structure 

5 . Overall gearbox dynamics and acoustics models 

6 . Comprehensive gearbox design criteria for reduced noise and vibration 
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